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We study a composite multimode light-two-level atom system in a cavity. We show that coupling 
of the two-level atom to multiple modes of the light destroys the Mott phase of the composite 
system thus making the system less useful platform for developing concepts in quantum information 
processing. 
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INTRODUCTION 

"dv 

The search for interesting and potentially useful quantum-mechanical phenomena on a mesoscopic scale in condensed 
matter and atomic physics is a challenging task. The Jaynes-Cummings model describes the interaction of a single, 
quasi-resonant optical cavity field with a two-level atom. The optical nonlinearities generated due to the coupling 
between the atom and the photons leads to an effective photon-photon repulsion [l|, Q . However it was shown that 
Oh- photon-photon repulsion degrades in the presence of many atoms 0, Later, Rebic et al. [B| showed that the 
nonlinear interaction afforded by placing a single two-level atom inside a cavity would suffice for realizing photon 
blockade. Greentree et al. Q showed that by adding photons to a two-dimensional array of coupled optical cavities 
each containing a single two-level atom in the photon-blockade regime, a long-lived, strongly interacting dressed 
states of excitations (coupled atom-photons) are formed which can undergo at zero temperature a characteristics 
Mott insulator to superfluid quantum phase transition. Devices based on photon blockade mechanism has a strong 
fS |' potential to be useful for quantum computation. In particular because of Mott's phase robustness, devices based on 
, l } ', this effect at non-zero temperature has been suggested @. Motivated by such interesting developments in photon 
blockade schemes, we study in this work, the prospects of utilizing quantum phases of light in a single atom driven by 
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THE BASIC MODEL 



multiple different modes of photons for possible quantum devices. In particular, we consider a two-dimensional array 
of photonic bandgap cavities. Each cavity contains a single two-level atom, quasi-resonant with the multiple cavity 
modes. In the classical limit, an array of coupled photonic bandgap cavities has been described for novel waveguide 
applications SMl, nanocavity lasers Q and in the quantum regime a two-cavity arrangement has been proposed 
as a Q-switch lllj. Recently,the dynamics of linear arrays of coupled cavities containing four-level atomic systems 
was studied by 111 ] . In the quantum optics context, the generalized Jaynes-Cummings model where the transition is 
mediated by multiple different modes of photons was first studied by Dantsker [12 1. 
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Eigenvalues and Eigenfunctions 

To create an atom-multi-photon system whose dynamics we want to study, we consider a two-dimensional array 
of photonic band gap cavities as discussed earlier [6(. Each cavity contains a single two- level atom, quasi-resonant 
with the multi-modes of the cavity. To motivate the search for Hubbard-model-type interactions within the present 
system, we first discuss the eigenvalues, second-order coherence and the probability of finding the atom in the excited 
state in the presence of multimodes. 

The system we consider here is an effective two- level atom with upper and lower states denoted by |1 > and |0 >, 
respectively. In the multi-photon processes, some intermediate states \i > , i = 2,3.., are involved, which are assumed 
to be coupled to |1 > and |0 > by dipole allowed transitions. Let ojq, oj\, and u>i, denote the corresponding frequency of 
the atomic energy level |0 >, |1 > and |t >, respectively. We also denote uj as the transition frequency between states 
|1 > and |0 >. The atom interacts with the n cavity fields with frequencies f2,;, i = 1, 2..n, where Qi + Q2 + ■■■Q n — u>. 
Also we assume that detuning between the atomic transition frequency and any one of the n modes is non-zero. 
Under these circumstances, the intermediate states can be adiabatically eliminated and the effective Hamiltonian of 
the two-level atom can be written in the rotating-wave approximation as 
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Here a] and fij are the creation and destruction operator of the n th cavity mode respectively, and a z , <r_ and 
a + are the usual pauli spin matrices. In the bare state basis |l,m,n 2 , ...n n > and |0,ni + l,n 2 + 1, ..n„ + 1 >, the 
eigenvalues and the eigenvectors are calculated as 



A ± = K(n 1 ,n 2 , ..n n ) ± Q(n ll n 2l ...n n ), 



where 
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Here A = lo — ^\ Wj. The corresponding eigenvectors are 



where, 



|±,m,n 2 ,..n„ >= A„ i;n2v . nn |0,m,n 2 ,..n„ > +-B nii „ 2 ,.. n Jl, m, n 2 , ..n„ >, 
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Fig. 1 shows a plot of the dimensionless eigenvalues for the 3 mode case. The photon energies have been subtracted 
for ease of comparison. The photon number in the first mode (m) is allowed to change while keeping the photon 
number in the other two modes (n 2 , n^) constant. The on-site photonic repulsion is evinced by the increasing energy 
separation with n\. The atom-photon coupling constants </j have been taken to same for all modes. 
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FIG. 1: Dimensionless eigenspectrum (E = A* — optical energy ) for a single atom interacting with 3-modes of a high Q-cavity 
as a function of A/g (dimensionless atom-cavity detuning), centered around zero (A* = Yli=i t^ifii). The eigenspectrum splits 
into two branches corresponding to the dressed states, |+, m, ri2, ■■n n > (upper branch) and |— , m, 712, ■■n n > (lower branch). 
The two branches anti-cross at A/g = 0, with the splitting increasing with increasing excitation number ni = 5, 10, 15 (in order 
of decreasing thickness of the plots) 

Population probability and Coherence function 

The results of the previous subsection leads to the following matrix representation of the time evolution operator 
for the given manifold M(ni,n2, ■■n n ). 
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The time evolution operator can now be utilized to calculate the density operator of the system at any time t with 
an arbitrary initial condition p(0) as 



p{t) = U{t)p{Q)U ] {t) 
The expectation value of any operator can then be obtained easily as 



(11) 



< 0(t) >= Tr[p(t)0(o)} (12) 
We assume that at t = the density operator can be decomposed into its atomic and field parts, i.e, 



P (0) =p A (0)®p F (0), 

where the atom is in the ith energy eigenstate \i > and the fields are in a general state as 



(13) 
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Futhcrmore the reduced atomic density operator p A {i) and the reduced field density operator p F (t) are introduced 
by taking the trace oip(t) over the field states and over the atomic states respectively. Now if we denote P nil n 2 ,~n„ (0) = 
\Rm,n 2 ,..n n | 2 as the initial photon distribution, then the probability of finding the atom in the excited state P\(t) for 
\i >= 10 > is calculated as 
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FIG. 2: Probability Pi(t) (left plot) of finding the atom in the excited state and coherence function (right plot) as a function 
of gt. The initial field states are prepared in 3-modes and 4-modes of coherent states with mean photon number fii = 20 in 
each mode. 
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The initial photon distribution is assumed to be coherent which implies that R ni ,n 2 ,..n n is written as 



exp 



R 



n,(* 



W2 



ni ,ri2 ,. -n-n 



vTL 



(16) 



Here rii is the mean photon number in the i th mode. The degree of n th order coherence with zero-time delay is 
defined as 
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whose magnitude controls the n-photon transition rate. 

If we assume the atom to be initially in the excited state then mean photon number of the ith mode and the 
expectation value of the correlation of photon numbers for the n modes that is related to the degree of second order 
coherence is 



< tii(t) >=Tr[p(t)hi], 
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< hih 2 --.n n >= Tr[p(t)hih 2 ...h r , 
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The time evolution of the probability of finding the atom in the excited state (P\(t)) and the n order coherence 
function (g n ) in the presence of 3-modes and 4-modes of the coherent state are plotted in Fig. 2. As seen from the 
figure the quantum revival and collapse of the Rabi oscillations are more compact for the 4-mode case than the 3-modc 
case. The more interesting observation is the fact that Pi is less for the 4-mode case than the 3-mode case. Similarly 
the coherence function g 4 < g 3 . This result has some interesting consequence for the quantum phase diagram of the 
composite photon-atom system to be calculated in the next section. 



QUANTUM PHASES OF LIGHT 



The Hamiltonian of our system is given by a combination of the Jaynes-Cummings Hamiltonian with photon 
hopping between cavities and the chemical potential term. Here we assume that only one mode (k th ) mode is able 
to hope between the cavities. The hopping is achieved by evanescent coupling between the cavities. The tunneling 
frequency is approximately given by tt — Qk/Q(£lk is the frequency of the k th mode and Q is the quality factor of 
the cavities ). If fi/. is kept large as compared to the frequencies of the other mode then one can achieve tunneling of 
only the k th mode. 



H = E H " c - K E 4 A ;±i - » E E + a i°i ( 2 °) 
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Here I is the site index. We are assuming nearest neighbour hopping and the hopping frequency k is assumed to be 
same. The chemical potential is same for all site if we assume zero disorder. Here ^li^id + a i~ a l\ IS * nc total 
number of atomic and photonic excitations (the conserved quantity in our system). We introduce a superfiuid order 
parameter ip =< a/ >, which we take to be real and use the decoupling approximation, a\a m —< a, > a m + < di > 
®m~ < ®i ><; >• The resulting effective mean- field Hamiltonian can be written as a sum over single sites 

h mf = ^2 H ? c - zK^^h + <i) + z ^\ 2 - E ( E + 4 -A (2i) 
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Here z = 3 is the number of nearest neighbours. To obtain the system's zero-temperature phase diagram, we use 
the procedure of ref. (l3j . All energies are now scaled with respect to ZK.The unperturbed ground state energy of the 

state with exactly a\ jo^j + crfaf^J particles is A = A~ — /! a\ + a^a^ , where A~ = \~ / zn. Wc only 

need to consider the negative branch for the purpose of determining the ground state since A~ < A + . A change in 
the total number of quasi-excitations per site will occur when A~ +1 — /2(n/c + 1) = A~ fc — fink- We can determine the 
critical chemical potential, p, c (rik), where the system will change from n k to n k + 1 quasi-excitations per site as 

Mc(n) = A- +1 - A- , (22) 
The second order correction to the energy with the well known expression 

E {2) = | < -, mi, m 2 , ..m n \V\-,n 1 ,n 2 , ..n n > | 2 ^ 23 ^ 
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Here we take the unperturbed eigenvector |— , nx,ri2, ■■n rl > corresponding to ^\ rij particles. Since the interaction 
V couples only to states with one more or less excitations than in the ground state, we find 

E (2) = (V^i£t + B- k B- k _^m k - l) 2 + (A~ k A; nk+1 ^m k + 1 + B mk B~ k+1 ^) 2 

\n k — \n k -l ~ P A m k — \n k + l + M 

According to the Landau procedure for second-order phase transitions, we write the ground state as an expansion 
in ip 



E m {ip) = a (m, p., A) + a 2 (m, fi, A)V> 2 + 0(ip 4 ), 



(26) 
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FIG. 3: Phase diagram of the coupled multimode light and the two-level atom syatem in a cavity. The left plot is for 3-mode 
case and the right plot is for 4-mode case. Dominating the extreme right-hand and the extreme left-hand edge (where photonic 
repulsion dominates at higher detuning) is the Mott-insulator (denoted by MI) and the superfiuid phase is found in between 
the MI regions (denoted by SF).The size of the Mott region is found to decrease with the quasi-excitations for both the 3-mode 
and 4-mode case. The size of the Mott lobes is found to be smaller for the 4-mode case as compared to the 3-mode case. 

E m (ip) is minimized as a function of the superfiuid order parameter tp. We find that ip = when 02 > and that 
ip =/= when a 2 < 0. This means that a 2 — signifies the boundary between the superfiuid and insulator phases of 
light. This yields 
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The subscript ± in Eqn. (27) denotes the upper and lower halves of the Mott insulating regions of phase space. 
Fig. 3 shows the plot of Eqn. (27) for the 3 mode case (left plot) and 4 mode case (right plot). By equating /2_ and 
/j+ we can find the point of largest A max for each Mott region (MI). The dynamics illustrated in Fig. 3 is extremely 
rich. Dominating the extreme right-hand and the extreme left-hand edge (where photonic repulsion dominates at 
higher detuning) is the Mott-insulator (denoted by MI) and the superfiuid phase is found in between the MI regions 
(denoted by SF).The size of the Mott region is found to decrease with the quasi-excitations for both the 3-mode and 
4-mode case. The transition from the SF to the MI phase of the quasi-excitations occurs at A max . The size of 
the Mott lobes is found to be smaller for the 4-mode case as compared to the 3-mode case. This observation can be 
correlated with the observations noted in Fig. 2. For the 4-mode case, the probability of transition to the excited state 
is less as compared to the 3-mode case. The photonic repulsion is a nonlinear process and is directly dependent on the 
probability of transition to the excited state. A lower transition probability implies a decreased photonic repulsion 
and hence a smaller Mott region. To find the cigcnenergies and to experimentally identify the various phases of the 
coupled atom-photon system, one can perform a transmission spectroscopy with the scattered light by direct read out 
of the number of photons coming out of the cavity. Photon loss can be minimized by using high Q cavities and thus 
ensuring that the light field remains quantum-mechanical for the duration of the experiment. 



CONCLUSIONS 



In summary, we have studied the strongly correlated dynamics of a two-level atom coupled to multi-mode cavity 
photons in the photon blockade regime. We have shown that coupling of the two-level atom to multiple modes of the 
light destroys the Mott phase of the composite system thus making the system less useful platform for developing 
concepts in quantum information processing. 
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